Order in a Spatially Anisotropic Triangular Antiferromagnet 



Scdigh Ghamari, 1 Catherine Kallin, 1 Sung-Sik Lee, 1,2 and Erik S. S0rensen 1 

1 Department of Physics & Astronomy, McMaster University 
1280 Mam St. W., Hamilton ON L8S J.M1, Canada. 

2 Perimeter Institute for Theoretical Physics 
31 Caroline St. N., Waterloo ON N2L 2Y5, Canada. 
(Dated: November 22, 2011) 

The phase diagram of the spin- 1/2 Heisenberg antiferromagnet on an anisotropic triangular lattice 
of weakly coupled chains, a model relevant to CS2CUCI4, is investigated using a renormalization 
group analysis which includes marginal couplings important for connecting to numerical studies of 
this model. In particular, the relative stability of incommensurate spiral spin density order and 
collinear antiferromagnetic order is studied. While incommensurate spiral order is found to exist 
over most of the phase diagram in the presence of a Dzyaloshinskii-Moriya (DM) interaction, at 
small interchain and extremely weak DM couplings, collinear antiferromagnetic order can survive. 
Our results imply that CS2CUCI4 is well within the part of the phase diagram where spiral order is 
stable. The implications of the renormalization group analysis for numerical studies, many of which 
have found spin-liquid like behavior, are discussed. 



I. INTRODUCTION 

Frustrated magnets are of considerable interest be- 
cause of their potential for exhibiting novel ground states 
in two or higher dimensions. 1 The spin- 1/2 Heisenberg 
antiferromagnet (HAF) on an anisotropic triangular lat- 
tice is a particularly well-studied model of frustrated 
magnetism and is relevant to CS2CUCI4, a material whose 
neutron spectrum has been taken as evidence for close 
proximity to a two-dimensional spin liquid phase. 2-4 Be- 
low 0.6K, CS2CUCI4 orders into an incommensurate spi- 
ral spin density state and is well described by the spin- 
1/2 HAF on a triangular lattice, with interchain diagonal 
exchange J' weaker than the intrachain exchange J, al- 
though weak intcrlayer and Dzyaloshinskii-Moriya (DM) 
interactions also play a role in stabilizing the ordered 
state. 

The spiral phase with ordering wave vector q c i = 
it + 2 sin -1 ( J'/2 J) is the classical ground state for the 
HAF on a triangular lattice for all J' < 2J. It is well 
established that the nearly isotropic (J 7 < J) spin-1/2 
HAF exhibits spiral order that smoothly connects to the 
three sublattice Neel order at J'=J. Several studies have 
proposed that, as J' is further reduced, quantum fluctu- 
ations destroy this order and stabilize a two-dimensional 
spin liquid phase. 4-10 Although this model is in proximity 
to a one-dimensional spin liquid (J'=0), it lacks the fea- 
tures typically associated with possible spin liquid order 
in higher than one dimension. That is, it neither exhibits 
macroscopic degeneracy of the classical ground state, as 
occurs for Kagome or pyrochlore lattices, nor substantial 
ring exchange, as occurs near a metal-insulator transi- 
tion. Nevertheless, essentially all calculations show that 
the quantum fluctuations are unusually large and the lo- 
cal moment has been estimated to go to zero for J'/ J as 
large as 0.9. 10 ~ 13 

While it has typically been assumed that the loss of 
spiral order signals the appearance of a quantum disor- 
dered state, Starykh and Balents, 14 using a renormal- 



ization group (RG) approach, found a collinear antifer- 
romagnetic (CAF) state stabilized at small J'/ J. 14,15 
They identified quantum fluctuations of order (J'/ J) 4 
that couple second nearest-neighbor (nn) chains antifcr- 
romagnetically, leading to CAF order, whereas a ferro- 
magnetic second nn chain interaction would be compat- 
ible with spiral order. This result is surprising since, 
although quantum fluctuations will typically select a spe- 
cific state from a degenerate manifold of classical ground 
states, it is unusual for fluctuations to stabilize an or- 
dered state of higher classical energy. 16 This order is 
rather subtle since it is only selected at 0(J'/J) 4 , so nu- 
merical studies taken as evidence for spin-liquid behavior 
may also be compatible with such order, as well as with 
spiral order which is strongly renormalized by quantum 
fluctuations. Exponentially weak spiral order has been 
predicted within a random phase approximation 17 which 
does not include the (J'/ J) 4 antiferromagnetic fluctua- 
tions identified in Rcf. 14. Numerical studies on multi- 
leg ladders find incommensurate spiral correlations with 
a wavelength which is strongly renormalized from the 
classical value and smaller than the finite system size for 
0.5 < J' I J < l. 18 This provides evidence for an incom- 
mensurate spiral ground state in this range of J', but one 
would need to study even larger systems for smaller J' 
since the spiral wavelength, A = 2n/(q — 7r), is expected 
to grow rapidly with decreasing J' /J. 

Calculations that directly compare the energies of the 
CAF and spiral states have given conflicting results. Par- 
dini and Singh, 19 using linked cluster series expansions, 
found that the spiral phase, with a variational order- 
ing wave vector, appears to have lower energy than the 
collinear state for J'/ J > 0.1. They also found the or- 
dered moment vanished rapidly with decreasing J'/ J, as 
expected for cither ordered state. By contrast, a cou- 
pled cluster method yielded a lower energy for the CAF 
state for J' / J < 0.6, although this method was unable 
to capture the vanishing moment. 20,21 Neither method is 
strictly variational, making direct comparisons of the en- 
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crgies difficult to interpret, as discussed in Ref. 19. Nev- 
ertheless, both methods show that the energy difference 
between these two states at small J'/J is significantly 
smaller than the energy scale set by quantum fluctua- 
tions, i.e., than the 0(J' 2 / J) difference between the clas- 
sical and quantum energies. 

Here, we reexamine the RG approach of Ref. 14, but 
rather than starting by integrating out every second 
chain as was done in Ref. 14, we adhere more strictly 
to the real-space RG method and keep quantum fluctu- 
ations of 0(J'/J) 2 . These are generated by marginal 
(and irrelevant) couplings and are not expected to or- 
der the system as J' — > 0. However, these fluctuations 
dominate the energy and are important for interpreting 
numerical results as they extend over an unusually large 
length scale (which diverges as J' — > 0). Surprisingly, 
these are of the opposite sign to the fluctuations respon- 
sible for CAF ordering. Thus, our RG analysis elucidates 
the competition between spiral and CAF order. We also 
present numerical studies that show that ferromagnetic 
fluctuations dominate at all J' < J, and, consistent with 
the RG analysis, these fluctuations only grow weakly with 
system size for small systems. 

The remainder of this paper is organized as follows. 
The model and RG equations is presented in Sec. II 
and the general RG flows and resulting ground states 
for varying initial conditions are discussed in Sec. III. 
Section IV discusses the specific initial conditions appro- 
priate for the Heisenberg model and the resulting RG 
analysis. The crossover from CAF to spiral order with 
increasing J' is also discussed in Sec. IV. Numerical stud- 
ies of this system are presented in Sec. V and compared 
to the RG results, while Sec. VI considers the effect of 
a weak Dzyaloshinskii-Moriya interaction. A summary 
and discussion of the results are given in Sec. VII. An 
approximate analytic solution to the RG equations, used 
in the analysis of the numerical studies, is given in Ap- 
pendix. 



where n indexes sites on the horizontal chains, y is the 

chain label, and J and J' are, respectively, intra- and 

interchain couplings as shown in Fig. la. The last term 

in Eq. (1) is the interchain DM interaction, where D = 
D2 .ii,i5 

We first consider the case D=0. The continuum ap- 
proximation is made for the horizontal chains, while the 
other direction is kept discrete. The spin operators are 
written in terms of the uniform and staggered magneti- 
zations: 22 S„ t y -> M- y (x) + (-l) n N y (x), where the sub- 
script y indexes the chains, and we use the sign conven- 
tion shown in Fig. la. 
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FIG. 1. (a) Triangular lattice with anisotropic exchanges J 
and J'. Staggered magnetization operators, N, at sites with 
the same color are accompanied with the same sign consis- 
tent with the depicted choice of y axis, (b) Phase diagram 
projected onto the plane of initial values of the relevant cou- 
plings, qm — 7*to- Dimer or other ordered states may be stable 
in the shaded region. 

The Hamiltonian for each chain is the Wess-Zumino- 
Novikov-Wittcn (WZNW) SU(2) model which defines 
the fixed point of decoupled chains, perturbed by intra- 
chain backscattering, with coupling 7b s : 23 



= + T*» I dx J R , y ■ J L J . (2) 



The interchain part of 7i in the continuum limit becomes: 



II. MODEL AND RG EQUATIONS 

Here, we derive the RG equations following an ap- 
proach similar to Ref. 14, but starting from short length 
scales and including marginal couplings and fluctuations 
of O(J'fJ) 2 in order to connect to finite size studies 
and to study the relative stability of collinear, spiral and 
dimer ordered states, as summarized in Fig. 1. Our start- 
ing point is the Hamiltonian: 



H — |^ S„ y ■ S n+ l y 



(1) 
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H into = j dx{^ M My ■ M y+1 + 7t 
y J 

(N, • d X N y+1 - d X Ny ■ Ny +1 ) +■■■} 



(3) 



where ellipses indicate irrelevant terms with more deriva- 
tives. M y is the sum of the left- and right-moving cur- 
rents, ^L,y and Ji?,y, which have scaling dimension 1, so 
My • My+i and backscattering both have scaling dimen- 
sion 2 and are marginal. The scaling dimension of N is 
1/2 and, therefore, that of the twist term (with coupling 
Itw) is 2. 

In addition to the above couplings, there are two 
relevant and two marginal interactions that are not 
prohibited by any of the symmetries of this model: 
<7ArN, y ■ N y+2 and g e e y e y+ 2 with scaling dimension 1, and 
•y £ (-l) y (syd x e y+1 -d x eySy +1 )/2 and g M M y -M y+2 with 
dimension 2. e is the staggered dimerization operator 
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and has scaling dimension 1/2. 24 Third and further nn 
chain couplings can be ignored for small J' / J because of 
their much smaller initial values. 

The theory is regularized by imposing a short distance 
cutoff for how close to each other the operators that per- 
turb % WZNW can be. Real-space RG along x is then con- 
veniently performed using operator product expansions 
(OPE) of the above interaction terms at separate space- 
time points and integrating over short relative spatial 
and temporal separations. 25,26 The OPEs of chiral cur- 
rents can be directly derived from the OPEs of chiral 
fermion fields, but the correct OPEs of N and e require 
bosonization. 22 Neglecting terms of 0(J' J J) 5 and higher, 
the /3-functions for the relevant and marginal coupling 
constants are: 27 

dilbs = lis ~ 6g 2 N , dilM = 7m ( 4 ) 

dlltw = -\lbsltw + iMltw - iltw9N - \lMle (5) 

dig N = 9n - \lbs9N + 9m(n + \ll w + 9m9n (6) 
di9e = 9e + ^lbs9e - f ffjuGv + \ll ~ %9m9n (7) 

dlje = f 7bs7e - ^ItwlM ~ \le9e (8) 
di9M = 9m - ^IbsChl (9) 
9iCn = -(n - |7fcs0v - Itw + 9m(n (10) 

diC M = ~2( M - 8^ 2 7 2 f + 76s Cm (11) 

All couplings are scaled by 2ttv = tt 2 J, and e l = L/Lq, 
so that l(L ) = 0, where L is the physical length and 
Lq is the initial length scale, which has been set equal 
to one in the above equations. Equations. (9)-(ll) are 
only brought in because these couplings affect the initial 
values of more relevant couplings, (jv is the coupling for 
the interaction term d x N y ■ 9 x Nj /+ 2, and C,m is that of 
(M y ■ M. y+2 ) (Jl, s +i • Jr :V+ i). Both of these are irrele- 
vant couplings with scaling dimensions 3 and 4. The f3- 
functions of jm and gu only describe the renormalization 
of the interchain backscattering part of M. y ■ Mj,/ (cou- 
pling between currents with opposite chiralities) which 
also are the parts that enter in the other /3-functions. 

The /3-functions of Ref. 14 differ from Eqs. (4)-(8) in 
that only the first three terms in the /3-functions for 
and g e [in Eqs. (6)-(7)] were included in Ref. 14 and 
the renormalization of marginal couplings other than lbs 
were not considered, itw, in particular, is needed to in- 
terpret the numerical results and to study the relative 
stability of the spiral and collinear AF states. As dis- 
cussed below, the 7 2 TO term in Eq. (6) is not expected to 
order the system, but describes the dominant fluctuations 



at short and intermediate distances and is important for 
interpreting numerical results on finite systems. The key 
RG equations for our study are Eqs. (4)-(6). 



III. GENERAL RG FLOWS 

At the lattice scale, only couplings present in the bare 
Hamiltonian, j tw = 7A//2 = J' /it 2 J. jbs = —0.23, 28 are 
nonzero. Integrating out short wavelength fluctuations 
up to the initial length scale L (larger than but compara- 
ble to the lattice spacing) generates all couplings allowed 
by symmetry, in particular, the two relevant couplings, 
g,/v(0) and <? e (0). Once generated, these couplings tend 
to grow exponentially in I (or linearly in L) while the 
marginal and irrelevant couplings grow at most logarith- 
mically with L. It follows from Eqs. (6) and (7) that, in 
the presence of negative backscattering, gM grows faster 
than g t . Therefore, for small J' /J, the value of <?jv(0) 
largely determines the fate of the system, as it will reach 
one (or J in unsealed units) at some length scale, I*, 
while the other couplings will remain much smaller. 

If 3at(0) is positive (i.e., antifcrromagnetic), it follows 
from Eq. (6), that <?at(/) remains positive and flows to +1 
at some I*. In this case, second nn chains will order anti- 
ferromagnetically. This is the result found in Ref. 14, 
where every second chain was integrated out to give 
g N (0) = -2g e (0)/3 = A$J' 4 /tt 6 J 4 , where 0.13 is 

a normalization factor. The CAF state is then stabilized 
for small J' / J, by the order-from-disorder mechanism, 
when one includes the effect of 7 m('*) ~ 0(J' / J). u 

However, numerics on finite systems show that second 
nn chains are coupled fcrromagnctically, not antiferro- 
magnetically as implied by gjv(0) > 0. Figure 2 shows 
the interchain Neel susceptibility, where a small finite 
staggered magnetic field is applied to one chain and the 
response of the second neighbor chain is calculated using 
exact diagonalization (ED) for three chains of length 12 
sites or less, with open boundary conditions. For a stag- 
gered field, h, along z, the response (Sf)/h, where site i 
is one of the center sites on the chain, is calculated. This 
response is the interchain Neel susceptibility, 



xs(l) = / dte(t)([si y (t),s* ty+2 (o)}) . 

71=1 J 

(12) 

It is found that (Si) aligns fcrromagnctically with the 
second neighbor chain for all < J' < 1, as well as 
for small J' < 0. Xs, which is discussed further in 
Sec. V below, is analytic in the couplings and, for small 
J' /J, is proportional to 3jv(L) to leading order. There- 
fore, it follows from these ED results that <7jv(0) = 
a(J'/J) 2 + b(J'/J) 3 + c(J'/J) 4 , where a, b, and c are 
all ferromagnetic and, in fact, ferromagnetic fluctuations 
dominate in these systems of up to 36 spins. \s calcu- 
lated for four chains with periodic boundary conditions 
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applied perpendicular to the chains (along y) yields ex- 
actly the same quadratic and cubic terms and a larger 
ferromagnetic quartic term. The difference in the quar- 
tic term is due to the longer range analog of jtw which 
couples third neighbor chains and is not included in the 
above RG equations as it has no qualitative effect on the 
running couplings. 27 
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FIG. 2. The response to a small staggered magnetic field, h, 
applied to one chain is studied using exact diagonalization for 
three chains of lengths from L = 6 to L = 12. The response 
of a central spin on the second-neighbor chain, {Si)/hL, is 
shown along with polynomial fits to the data. The sign of 
the response is found to be such that (Sf) always aligns fer- 
romagnetically with the other chain. 



The calculated \s is sensitive to boundary conditions 
for small systems. Periodic boundary conditions (PBC) 
along the chains (along x) frustrate spiral correlations 
generated by Jtwi which we will see in the next section 
are connected to ferromagnetism. Indeed, for the small 
systems studied with PBC along x (up to 30 spins, not 
shown here) we found Xs is of order (J'/ J) 4 and has the 
opposite sign to the data in Fig. 2. However, the fact that 
spins Si tV and 5j_|_ XiJ ,-|_2 for small s are coupled ferromag- 
netically at order (J'/ J) 2 and the higher one is indepen- 
dent of the boundary condition and is seen even in small 
systems with PBC along x. In particular, (Si jV ■ Si t y+2) 
is always positive irrespective of the boundary condition. 
This implies that gjv(0) < 0. In fact, for open boundary 
conditions, our ED studies find that (— l) x {Si V -Si+ x ,y+2) 
is positive for all x, all systems sizes and J' values studied 
(J'/J < 0.5). 

It follows from the RG equations that for <7at(0) < 0, 
9n{1*) can be cither negative or positive, depending on 
the precise value of 3at(0). In this case, the flow of gm is 
sensitive to the magnitude of the initial value due to the 
competition between <?at(0) < and terms in Eq. (6) that 
drive <7jv(0 antiferromagnetic. This sensitivity to initial 
conditions can be seen by numerically studying the RG 
flows for an initial <7at(0) = a^f w (0) as the constant a 
is varied. Here, we use bare values (i.e., values at the 



lattice scale) for the initial values of the other couplings. 
RG flows are shown in Fig. 3 for -0.30 < a < -0.26. 
For a > —0.28, one finds that gN flows to 1 (strong an- 
tiferromagnetic coupling) and jtw remains small, while 
for a < —0.28, gN flows to -1 (strong ferromagnetic cou- 
pling) and jtw increases marginally. As discussed below, 
g n — > — 1 signals spiral order, whereas gN — > 1 signals 
CAF order as studied by Starykh and Balcnts. 14 How- 
ever, as seen in the inset of Fig. 3 for gN — > 1, even 
in the case of CAF order, the initial RG flows can dis- 
play increasing ferromagnetic coupling between second 
nn chains. In Sec. IV, it will be shown that this is the 
case for the Heisenberg model. In a narrow range for 
a —0.28 (which widens as J' is increased), the growth 
of gN is hindered and j tw , or g e reach unity first. This 
crossover value, g N lt , depends on 7h s and J'. Depend- 
ing on initial conditions, columnar or staggered dimer- 
izations or more complicated incommensurate states can 
be stabilized in this crossover region, as denoted by the 
shaded area in Fig. lb. 

In the ferromagnetic regime where gN grows exponen- 
tially to negative values (g N < g™'), sites along the y 
direction align, consistent with zero y component of the 
ordering wave vector. In this regime, jtw, which stabi- 
lizes spiral order, grows marginally. In the classical limit, 
the ordering wave vector q = (it + e)x is (ir + J'/ J)x for 
J' « J. The effect of quantum fluctuations on q, or 
e, can be determined from the RG analysis. The most 
robust spiral state occurs if ferromagnetism is selected 
at quadratic order in J'. The length scale, I*, at which 
gN becomes comparable to J in magnitude is inversely 




Ytw 

FIG. 3. Flows of (/jv(0 and ytm(l) for several initial values 
with ferromagnetic <?iv(0) = c*7t™(0), 5e = 7e = 0, and J' = 
0.05. gN and jt w are in scaled units (divided by it 2 J), and 
a is varied from -0.26 (CAF flow closest to vertical axis) to 
-0.30 (spiral flow closest to vertical axis). The inset zooms in 
on the shaded rectangle. Flows emanating from initial values 
very close to g^ 1 * have flow lines substantially away from the 
ftw = axes and grow considerably slower than those close 
to the vertical axes. Note the extreme sensitivity to initial 
values close to the critical value g^'* as shown in the inset. 
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proportional to the initial value of g^ which, in turn, 
is oc J' 2 in this case. Since 7^ changes logarithmically 
and remains negative as g^ grows, at this length scale, 
we expect almost the same intrachain Hamiltonian but 
with a somewhat renormalized j tw ■ This Hamiltonian de- 
scribes interactions between blocks of Neel-ordered spins 
of length I* and can be treated classically. Thus, the 
ground state has a spiral order with e oc (J'/ J) 3 . Since 
1/e ^> I* for small J', the pitch of the spiral is much 
longer than the Neel blocks and this treatment is justi- 
fied. This is a significant renormalization of the ordering 
wave vector toward the ID limit of 7r, where spins on 
neighboring antiferromagnetic chains are oriented at tt/2 
with respect to each other and requires that ferromag- 
netism is selected at quadratic order. 

Even more fragile spiral order is stabilized if Ag = 
9n(Q) — ff™' ~ an d ferromagnetism is selected at cubic 
or quartic order in J'. Then e oc (J'/J) n , where n = 4 
(or 5) for cubic (or quartic) selection. The exponentially 
weak spiral order, e ~ e ~ a ( J / J ) found within a random 
phase approximation, 17 also follows from the RG equa- 
tions if one assumes that neither g^ nor g e can grow. 
To leading order in J'/ J, Eq. (5) describes this expo- 
nentially weak spiral order, where a = 2.6 if the value 
lbs ~ —0.38 at ao is used. However, since there is no 
symmetry that prevents gw (or g e ) from growing at all 
orders, one does not expect an exponentially weak spiral 
order to be stable for very small J'/ J. 

As J' increases, even for initial conditions that favor 
antifcrromagnetism, Ag = 5at(0) — c/™* < 0, the RG flows 
move toward increasing j t w, i-e., toward the spiral state. 
This is because jtw is boosted by both backscattering 
and YjVf- The latter enhances the spiral incommensura- 
tion, e, so it is larger than e ~ a ( J / J ) . For the Heisenberg 
model, where 7b s (0) is sufficiently negative to suppress 
the growth of g t relative to gN, the competition is pre- 
dominantly between the CAF and the spiral states. If a 
sufficiently strong frustrating second nn interaction along 
the chains is added to the Hamiltonian, the competition 
is then between the CAF and dimerized phases as dis- 
cussed in Rcf. 14. 



IV. INITIAL CONDITIONS AND RG RESULTS 

It follows from the above discussion of the RG flows 
that it is important to determine <?at(0) accurately. Nu- 
merics on small systems clearly show that <7at(0) is fer- 
romagnetic, whereas Eq. (6) would naively predict that 
9n{0) ~ +7t 2 u;- I n fact, only antifcrromagnetism appears 
to enter the RG equation for g n to all orders in J'/ J. 
This would imply not only that spins on second neigh- 
bor chains are always antiferromagnetically ordered, but 
that the dominant fluctuations at all length scales are an- 
tiferromagnetic. Actually, one sees how second neighbor 
chains can be coupled ferromagnctically by considering 
the effect of a twist in one chain on spins in neighboring 
chains. Treating 7tu,Nj,-|-i • d x lSSy as perturbations to the 



decoupled ID chains, one finds in second order that fluc- 
tuations in the y-th chain generate a local ferromagnetic 
coupling, H' a ~ -i 2 w (\d x Ny\ 2 )N y ^i ■ N y+1 . Although 
the static twist (c^Ny) is small, the fluctuation can be 
O(l), mediating an interaction of order J' 2 . 

The peculiar behavior where the sign of the quantum 
correction for g^ depends strongly on scale can also be 
seen in the continuum theory, once the theory is regu- 
larized. The OPE between two twist terms generates a 
term in the effective action: 



- 7 2 w j dzdz , [d x d x ,G(z-z , )]Ny- 1 (z)>-Ny +1 (z'). (13) 

Here z = i(vt - x), G(z - z') = (N y (z) ■ N„(z')) ~ 
(z — z') 2 + a 2 , and a c is larger than but compara- 
ble to the lattice spacing. The integrand in Eq. (13) is 
positive for < \z — z'\ < a c , which generates the ferro- 
magnetic initial condition for g^, but becomes negative 
for \z — z'\ > a c , favoring antiferromagnetic correlation 
at longer lengths. 

The presence of a ferromagnetic coupling at short dis- 
tances is very general and independent of the specific reg- 
ularization scheme: it simply follows from the fact that 
[d x d x >G(z — z')] in Eq. (13) is a total derivative which 
integrates to zero and is negative at large distances. The 
presence of this total derivative implies that this term 
does not contribute at q = in a gradient expansion and, 
consequently, is not expected to order the system. 29 For 
example, if one fully integrates out every second chain to 
generate effective couplings between second nn chains, 
as was done in Ref. 14, this term vanishes. This implies 
that the initial value of <?at at quadratic order must be 
such as to cancel the effect of j 2 w in the /3-function in 
the q — > limit. This places constraints on the initial 
conditions that are discussed below. In general, how- 
ever, <?at(0) = aJ' 2 + bJ' 3 + cJ' 4 , where a, b, and c are 
all negative (i.e., ferromagnetic). The cubic and quartic 
contributions to the initial condition are related to the 
quadratic and cubic contributions in the beta function 
for jtw 

Wc note here that the jtw term in the /3-function for g^ 
is not the only term that contributes ferromagnetism at 
short length scales. For example, the longer range analog 
of itw and irrelevant couplings such as (at (or even higher 
derivative terms) contribute to the /3-function for g^ in a 
very similar way to the j 2 w term, but with smaller coef- 
ficients. While the inclusion of such irrelevant terms can 
have a quantitative effect at short lengths, the qualita- 
tive behavior of strong ferromagnetism, which persists to 
fairly long-length scales, is captured by including the ef- 
fects of the marginal coupling, jtwj only and is unchanged 
by the inclusion of more irrelevant couplings. 

That gN should not grow exponentially due to the jf w 
term is not only suggested by Eq. (13), but can be seen 
by expanding the interchain susceptibility, Xs, hi pow- 
ers of J', treating the interchain interaction as a per- 
turbation. One finds, using the simple power counting 
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suggested by scaling, that Xs{L)/L grows at most loga- 
rithmically in L at quadratic and cubic orders, but there 
is a contribution at quartic order that grows linearly in 
L. It follows from this perturbative expansion that the 
jf w term causes no exponential growth in <?at(Z)- The 
initial condition that corresponds to "tuning" (no expo- 
nential growth) is given in Eq. (A. 9) in Appendix, where 
the RG equations are solved analytically to cubic order. 
For 7b s = and the initial condition, </at(0) = — J' 2 /4, it 
follows from Eqs. (5) and (6) that <7at(7) remains constant 
at quadratic order. However, for nonzero 7& s , 5at(0 flows 
even at quadratic order and a larger ferromagnetic ini- 
tial condition is needed to ensure that g^ does not grow 
exponentially in / at quadratic order. 

While there is a ferromagnetic quartic contribution to 
<7at(0) from the j% w term, there is also an antiferromag- 
netic contribution from the gmCN term in the /3-function. 
This term, generated within the operator product ex- 
pansion technique, is analogous to the antiferromagnetic 
initial condition identified by Starykh and Balents. 14 
gM^N is zero at the lattice scale since these couplings 
do not appear in the microscopic Hamiltonian, but is 
nonzero in subsequent RG steps and provides an ini- 
tial condition equal to 8-/ bs J' 4 /n 9 J 4 = 0.06 J' 4 /tt 6 J 4 at 
I = if one sets 7m (0) = 2ltw(0) = 2J'/ir 2 J. This 
differs from the Starykh and Balents initial condition 
9n(0) = AqJ' a /tt 6 J 4 by about a factor of two, but both 
calculations rely on the continuum approximation which 
is only approximate near the lattice scale. Also, the two 
results are almost identical if one uses the value of 7& s 
extrapolated to the lattice scale. An analysis of both 
calculations shows that they rely on the same physics, 
which involves generation of a gradient coupling between 
second- neighbor chains ((n) and requires chiral symme- 
try breaking (i.e., nonzero "fbs in the continuum RG lan- 
guage). 

The long-distance RG results of Ref. 14 are recovered 
within our analysis for sufficiently small J' if the ferro- 
magnetic initial condition associated with jf w is equal to 
g^ 11 to C(J'/J) 4 . In this case, the quartic term, <7mCat, 
provides the essential contribution to the initial condi- 
tion of gM that drives the ordering. For this initial con- 
dition, the RG flow for small J' corresponds to one of 
the flows shown in Fig. 2 where g m is initially ferromag- 
netic, but passes through zero at an intermediate length 
and then becomes large and antiferromagnetic. For small 
J', the intermediate length where g^ crosses over from 
ferromagnetic to antiferromagnetic, is proportional to 
(J/ J') 2 , neglecting logarithmic corrections. This follows 
from the fact that the quadratic contribution to <?./v(0) 
grows at most logarithmically in L, while the quartic 
antiferromagnetic contribution grows linearly in L. Con- 
sequently, the second nn chains in finite systems with 
chains of length L will be fcrromagnctically correlated 
for L < Lfm = A(J/ J') 2 . The coefficient A is estimated 
to be 20-40 lattice spacings, depending on which value is 
used for the antiferromagnetic contribution to <7jv(0)- 

A direct transition from CAF to spiral order is ex- 



pected to be discontinuous, while the RG flows change 
continuously. Nevertheless, one can extract some infor- 
mation about this transition from the RG analysis. As 
J' increases, the RG flows move toward larger y tw at 
/*, which shows a tendency of moving toward spiral or- 
der. From the numerical solution of the /3-functions of 
Eqs. (4)-(8), it follows that J' c ~ 0.3, where J' c is de- 
fined as the value at which ^ftwi}*) first reaches J while 
\9n{1*)\ < J- For 7 bs (0) = -0.23, g e (l*) remains small 
compared to J at this J', so the RG predicts that there is 
no intermediate dimer phase between the CAF and spi- 
ral ordered phases. This estimate for J' c results from us- 
ing the value A$ J' A /-k®J a for the AF contribution to the 
initial condition 14 and setting the ferromagnetic initial 
condition or tuning as discussed at the end of Appendix. 

Terms higher order than (J'/ J) and irrelevant cou- 
plings which are ignored, as well as exactly where one 
terminates the RG flows, can all affect J' c so J' c ~ 0.3 is 
only a crude estimate. The main reason for the somewhat 
small value of J' c is that the CAF state is only selected 
at quartic order, while the competing marginal coupling, 
jtw is linear in J' and boosted both by backscattering 
and by 7m, which is also linear in J'. We note that the 
estimate of J' c is compatible with the observation of spiral 
correlations at J' = 0.5 in numerical studies of multi-leg 
ladder systems. 18 

V. NUMERICAL RESULTS 

The ED results shown in Fig. 2 show only ferromag- 
netism at all J' , but the system sizes are not large enough 
to determine whether the strength of the ferromagnetism 
is changing with increasing system size for small J'. To 
address this, we turn to a more detailed analysis of the 
numerical studies of the interchain susceptibility, Xs, us- 
ing density matrix renormalization (DMRG) to study 
larger systems. Specifically, we address the question of 
what initial conditions are compatible with the study of 
finite systems or, alternatively, are finite size studies com- 
patible with the RG analysis. Using finite-size scaling, 
the susceptibility Xs(L,{gi(0)},L o ) of the system with 
size L = L$e l with couplings {<?i(0)} defined at scale Lq 
(I = 0) can be related to the susceptibility of the system 
with size Lq as 

X.(£,{ffi(0)},£o) = ^[l- lbs {L )l)] 1/2 Xs{L Q , { 9l (l)},L ), 

where the linear-L dependence in the prcfactor is due 
to Xs having scaling dimension 1 in the absence of the 
backscattering, and the second factor is the contribu- 
tion of the running backscattering that modifies the 
scaling dimension of Xs- For small <?jv(0) (small J'), 
Xs(Lq, {gi{l)}, Lq) is analytic in the couplings and pro- 
portional to <7at(0 to leading order, and we have Xs(L) oc 
(L/L )[l — 7bs(^o)0] 1 ^ 2 .9A r (0- There is no further quan- 
tum correction to this relation to leading order in J' /J, 
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since <?jv(0 ^ s fully renormalizcd and includes all quantum 
fluctuations. 

Since three chains captures all contributions to <?at and 
Xs to order (J'/ J) 3 as well as the key (J' /J) 4 fluctuations 
that drive CAF order, we confine our studies to three- 
chain systems. The interchain susceptibility is studied 
using ED and DMRG for the HAF on three chains with 
open boundary conditions for 24 to 84 spins (L=8 to 28); 
ED was used for up to 36 spins. Up to m = 2600, states 
were kept in the DMRG, while typically keeping the pre- 
cision e < 10 -8 . The results arc shown in Fig. 4. For all 
values of < J' < J, \s was found to be ferromagnetic. 
As noted earlier, ED studies also found that the corre- 
lation between second nn chains has quadratic, cubic, 
quartic (and higher order) in J' contributions which are 
all ferromagnetic for L < 12 (as expected from analyzing 
the term), and the additional quartic contribution 
from 4 chains is also ferromagnetic. There is no sign in 
the data of a turnover to negative values of Xs expected 
for AF correlations. As J 1 / J increases, the length scale 
at which antiferromagnetism should be observable if the 
CAF state is stable becomes shorter. For example, at 
J'/ J = 0.5, the RG analysis predicts that one should see 
a turnover in \s for L > 20, although we don't expect RG 
to be quantitatively accurate at such large J'/ J. As J'/ J 
increases, higher-order fluctuations that are not included 
in the RG and that are not fully captured by studying 
three chains, can either further stabilize or destabilize 
CAF order. 
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FIG. 4. The staggered Neel susceptibility, Xs> f° r three chains 
as a function of chain length, L, at different J' values, calcu- 
lated by ED and DMRG. The sign of Xs is such that second nn 
chains are ferromagnetically coupled for all L and J' studied. 

For L > 10 and J'\ < 0.1J, the change in Xs as a func- 
tion of L is fit to the predicted RG flow for g^ , Jtw , and 
jbs as functions of initial conditions, treated as fitting pa- 
rameters. 7m is incorporated approximately to capture 
the cubic behavior of <7jv(0 as described in the Appendix. 
The restriction to small J' is to ensure that Xs is propor- 
tional to gN and higher order terms can be neglected. 
Figure 5 shows Xs scaled by (L/L )[l - jbs{L )l)] 1/2 ! , 
with the lines being the RG fits to the data. One can see 



the different slopes of curves for positive and negative 
values of J', which shows the effect of the cubic term. It 
is also clear from these curves that g^ is growing only 
weakly with L. 

The initial conditions at L = 10 extracted from the 
numerical data are jtw = 0.416J' + 0.121J' 2 , g^ = 
-0.0450 J' 2 -0.0425 J' 3 , 7 hs = -0.072, and lM = 0.24 J'. 
While the extracted value of jbs(Lo) is substantially 
closer to zero than the value of —0.19 for periodic bound- 
ary conditions, this is consistent with an independent 
calculation of Xs for a single chain. For short chains, 
open boundary conditions are known to favor dimeriza- 
tion relative to periodic boundary conditions. 30 The fact 
that the ratio of 7(u,(0)/7m(0) extracted from the numer- 
ical data is noticeably larger than the ratio at the lattice 
scale can be attributed to the effect of irrelevant cou- 
plings that contribute ferromagnetism in an analogous 
way to the jf w term. 

Although the extracted initial conditions are such that 
<7jv flows to a large ferromagnetic value, which would 
seem to indicate spiral order, the deviation from tuning 
at both quadratic and cubic orders is sufficiently small 
to be compatible with finite size effects. For the values 
of 7ts(0), 7tu;(0), and 7m (0) extracted from the data, 
tuning corresponds to g N (0) = -0.0447 J' 2 - 0.0479 J' 3 + 
0{ J' / J) 4 . Therefore the quadratic (cubic) term deviates 
by 1% (12%) from the value predicted from the RG equa- 
tions for an order selected at C(J'/J) 4 . Given the un- 
certainty introduced by finite-size effects for 10 < L < 30 
and the approximations made in the RG, the numerics 
for small J' are consistent with no order being selected 
to cubic order, while the lack of any tendency toward an- 
tiferromagnetism with increasing L for J'/ J > 0.5 sug- 
gests the CAF state is not stable for these larger values 
of J'. The numerics also confirm the RG prediction that, 
not only is there ferromagnetism at the lattice scale for 
all values of J' / J, but this ferromagnetism continues to 
grow to fairly long length scales. Ideally, one would like 
to study the quartic term for small J' as a function of sys- 
tem size since the RG analysis predicts that CAF order is 
selected by an initial condition for g^ that deviates from 
tuning by about 20% at quartic order. However, this 
would require studying larger system sizes while keeping 
more states in the DMRG calculations and is beyond the 
scope of this study. 



VI. DZYALOSHINSKII-MORIYA 
INTERACTION 



It is also interesting to consider the effect of a nonzero 
DM interaction in this RG picture, since it is known that 
the interchain DM interaction further stabilizes spiral or- 
der in this model. 11,31 In the continuum limit, the inter- 
chain DM interaction can be written as e ab NyNy +1 with 
the initial condition g D (0) = 2||Dj| = 2D/tt 2 J. Note 
that this term breaks the SU(2) symmetry of the original 
model and will order the spins in a plane perpendicular 
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FIG. 5. The scaled susceptibility, \s = Xs/[L{1 -Jbs (0)l) 1/2 \, 
which is proportional to the relevant coupling, gM, for small 
J', is shown for different values of J' /J. Points represent data 
from ED (for L < 12) and DMRG (L > 12) and lines are fits 
to the RG equations, where the initial conditions are used as 
fitting parameters. 



to the direction of D. Just like the term, the 

DM interaction has scaling dimension 1 and is relevant. 
It flows according to the following /3-function: 



the absence of any DM interaction, the estimated DM 
interaction for CS2CUCI4, D w 0.05J, is much greater 
than what is required to stabilize spiral order. 

One can also consider the effect of an intrachain DM in- 
teraction with the Hamiltonian ^ n D' • (S n . y x Sn+i^). 
This interaction can promote several different phases 
even in a single spin- 1/2 chain for interchain DM cou- 
pling, D', comparable to J. 32 For D' <^ J, it has little 
effect on the phase diagram shown in Fig 6. One can show 
that the primary effect of D' is to renormalize the inter- 
chain DM coupling, D D-D'J'/2J. 33 For D = 0, this 
implies that D' > 0(J' 3 /J 2 ) will stabilize spiral order at 
small J'/ J. For nonzero D and |D'| ~ |D| or smaller, 
the intrachain DM coupling can slightly tilt the ordering 
plane and/or slightly shift the phase boundary between 
CAF and spiral, depending on the relative directions of 
D and D'. This can be understood from the fact that, 
for small J' /J, the CAF and spiral states exhibit essen- 
tially the same AF order along a chain and only differ 
significantly in how the chains order with respect to each 
other. 




digo = 9t> - -^IbsgD - 4g N g D + -7mPd ■ (14) 

The DM interaction also introduces a — 2g 2 ) term (for 
components perpendicular to D) in the /3-function of 
gM, Eq. (6). In addition to getting some boost from 
the marginal backscattering and 7m, the DM interac- 
tion both is enhanced by and promotes ferromagnetic 
gM- Consequently, it suppresses CAF order relative to 
spiral order. In any region of the phase diagram where 
spiral order is stable for D = 0, the DM interaction will 
further stabilize this phase relative to CAF order. 

Since both g^ and gjq have the same exponential 
growth for small J', a DM interaction with <7r>(0) greater 
than the untuned (antiferromagnetic) part of gjv(0) will 
overcome the CAF order as gu will grow to be of 0(1) 
first. The DM interaction favors neighboring Neel chains 
to be aligned perpendicular to each other, but in the pres- 
ence of jtw (which grows marginally), a spiral state will 
be stabilized. This suggests the phase diagram shown 
in Fig. 6 with a transition between the CAF and spiral 
states at D oc J' 4 / J 3 for small J'/ J. This transition is 
expected to be first order, so the RG crossover between 
these two states is only approximate. However, because 
CAF order is only weakly favored, at order (J'/ J) 4 , the 
maximum D for which CAF order survives will also be 
very small and is estimated from the RG analysis to be 
less than 10~ 4 J if J' r < 0.3 J. Even if the interchain cou- 
pling for CS2CUCI4 is such as to stabilize CAF order in 



FIG. 6. Suggested phase diagram where the critical inter- 
chain Dzyaloshinskii-Moriya coupling, D, vanishes as J' 4 /J 3 
at small interchain coupling J'. The ordering wave vector of 
the spiral state varies continuously and is q — (n, 0) for J' = 
and nonzero D. The region of stability of the CAF state is 
very small since 7?Max is estimated to be less than 10~ 4 J. 



VII. CONCLUSIONS AND DISCUSSION 

In summary, the full RG equations for the HAF on an 
anisotropic triangular lattice describe the intense compe- 
tition between CAF and incommensurate spiral orders. 
A tiny perturbation of C(J'/J) 4 can tilt the balance be- 
tween these two states at small J'/ J, while a larger per- 
turbation is needed to stabilize dimer order. As predicted 
by Starykh and Balents, 14 the CAF state is stable at suf- 
ficiently small J'/ J, but either an interchain DM interac- 
tion or a ferromagnetic second neighbor interchain inter- 
action of C(J'/J) 4 would stabilize an incommensurate 
spiral state. By contrast, an 0{ J) intrachain interaction 
that pushes backscattering, 7& s , toward zero or an inter- 
chain dimer coupling, which is larger than 0(J' /J) 4 is 
needed to stabilize dimer order. The fact that such an 
extremely weak interchain DM coupling is sufficient to 
stabilize spiral order at all J' /J, implies that CS2CUCI4, 
with an estimated D ps 0.05J, is well within the region 
of stability for spiral order and is not close to a CAF 
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instability. 
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as at J'/J ~ 0.9. 6,12,13 However, numerics which are in- 
terpreted as loss of spiral order could be compatible with 
the weak spiral order expected in the presence of quan- 
tum fluctuations. Furthermore, since the wavelength of 
the spiral grows rapidly with decreasing J'/J and its en- 
ergy lies so close to that of the CAF, as well as to many 
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L F m = A(J/J') 2 . The ED and DMRG results are consis- 
tent with this and suggest that the antiferromagnetism 
is hidden by larger ferromagnetic fluctuations except at 
very long length scales. One might be able to see these 
ferromagnetic correlations in the CAF state using cluster 
or series expansion calculations for the infinite system. 

While our numerical studies are consistent with the RG 
analysis, they do not provide any direct evidence of CAF 
order. Given the long- length scales involved for J' < J' c , 
it would be extremely challenging to identify this order 
from finite-system-size studies. For small J'/J, a large 
static magnetic susceptibility has been observed at or 
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7T, 
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but this feature is common to 



both CAF and spiral fluctuations. The two states would 
be distinguished by the q y dependence of the suscepti- 
bility and very long chain lengths are required for these 
correlations to develop. The interchain Neel susceptibil- 
ity, \s, together with finite-size scaling may provide one 
fruitful avenue for such studies. For example, one could 
simply look for any indication of a sign change in \s with 
increasing L for J'/J < 0.5. Our numerical studies are 
restricted to L < 30, because of the high precision re- 
quired to extract RG parameters and give no evidence of 
a sign change. However, to within the expected accuracy, 
the numerics are consistent with the RG prediction that 
quadratic and cubic (in J'/J) fluctuation effects do not 
select any order. 
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Appendix: Approximate Analytic Solutions 

To better understand the flow of g^ and j t w under the 
RG, one may attempt solving a reduced subset of the 
given /3-functions [Eqs. (4-8)] analytically. If one drops 
terms in the /3-functions that contribute at quartic and 
higher orders (in J'/J) to gx, partial analytic solutions 
can be obtained that yield insight into the initial condi- 
tions corresponding to "tuning", i.e., to no exponential 
growth in gjy. This leaves us with the following reduced 
/3-functions: 



dabs = ll s , da M = lli 



ilbsltw + iMlt w 



digN = gN - \ib S gN + \itw > 

which lead to the analytic expressions: 



76s (0 



76,(0) 
1-76 S (0)Z 



, 7m(0 



7m(0) 



I- 1m {0)1 



7ttu(0 = 7t<u(°) 
9n(1) 



1-7m(0)Z 

7L(0) f l e-yi- 76,(0)* 
4 Jo (l- 7M (0)t) 2 



(A.l) 
(A.2) 
(A.3) 

(A.4) 
(A.5) 



dt + g N (0) 



Vl- 76,(0)/ 



(A.6) 



For a non-zero jm (0), the integral in Eq. (A.6) cannot 
be expressed in terms of elementary functions. However, 
assuming 7 m(0)£ <C 1, it is possible to expand ,. - 1 
and find an expression for <7_/v(/) as series of lower, incom- 
plete gamma-functions: 



9n{1) 



^P-e WO) ^/3^(6) ]>> + 1)7^(0 

n 

-9n(0) eVl- 76,(0)/ (A.7) 



where F n (l) is given recursively by: 
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F n (l) = 



E 



,1- 



^y) -T L {n+ |,-^oy) 



(-1) 



^ 7 -(0)(n-m) 



. x n — rn 



(I) 



(A. 



Each term in Eq. (A. 7) is accompanied by a power of 
7i/(0), and also we have a 7^(0) factor that multiplies 
the sum. Since the above analytic expressions for the 
couplings were derived by neglecting quartic and higher- 
order contributions, we dismiss all terms in the expan- 
sion except the first two terms. Finally, using the above 
analytic expressions, the tuned initial condition for the 
Neel-Neel coupling, g^ % (0), is found [ 76s (0) = -0.23]: 



gTW = - 0.27620(J'/J) 2 



-1.20152(J7J) 3 +e>(j'/J) 4 
(A.9) 



It is also possible to determine the coefficient of the quar- 
tic term in the above expansion numerically, which gives 
8.0. Note that we have taken ^'(0) as the initial value 
for which gAr(7 final ) = [7tto(/ f i„ a i) = 1] which for any 
J' < 0.2 is not any different from other ways of defining 
tuning (e.g. givOfmai) = — !)• However, the estimated 
value, J' c / J , for the transition from CAF to spiral, would 
be smaller if <7jv(i final ) < was used. 



Finally, note that in the above calculation 76s (0) = 
—0.23 was used, which is the coupling of backscattering 
term at 4ao rather than ao where jm and jtw have their 
bare values. Of course nothing is special about ao and 
one can repeat the same analysis at any arbitrary length, 
knowing the values of the couplings at that length scale. 
Thus what has been neglected is the small growth in 
7m and jtw from ao to 4ao- Taking into account this 
growth of these couplings, the quadratic and cubic parts 
of 5^(0) become, respectively, 0.291374 and 1.31288. 
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